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We present simple lattice realizations of symmetry-protected topological (SPT) phases with g-form 
global symmetries where charged excitations have q spatial dimensions. Specifically, we construct 
d space-dimensional models supported on a (d -|- l)-colorable graph by using a family of unitary 
phase gates, known as multi-qubit control-^ gates in quantum information community. In our 
construction, charged excitations of different dimensionality may coexist and form a short-range 
entangled state which is protected by symmetry operators of different dimensionality. Non-triviality 
of proposed models, in a sense of quantum circuit complexity, is confirmed by studying protected 
boundary modes, gauged models and corresponding gapped domain walls. We also comment on 
applications of our construction to quantum error-correcting codes, and discuss corresponding fault- 
tolerant logical gates. 


I. INTRODUCTION 

The study of symmetry-protected topological (SPT) phases has attracted a considerable amount of attention [I]- 
[21] . Recently, generalizations of SPT phases with higher-form symmetry have been discussed [22U26] . Ordinary SPT 
phases are discussed in the presence of a global 0-form symmetry operator of the on-site form: 

C7(9) = (g) (1) 

3 

where g G G is an element of the symmetry group G and j represents a lattice site. The symmetry operator imposes 
a conservation law where charged excitations are point-like objects. A g-form global symmetry can be imposed by an 
operator of the form which acts on a closed codimension-g manifold M. (codimension g-h 1 for a space-time). 

In such a theory, charged excitations have g space dimensions and symmetry operators impose conservation laws 
on higher-dimensional charged objects. There have been several pioneering works in this direction [23H26) . Namely 
the work by Kapustin and Thorngren proposes a family of lattice realizations by replacing a finite group by a finite 
2-group 1^ . 

The goal of this paper is to present simple lattice realizations of bosonic SPT phases with higher-form global 
symmetry and discuss their quantum information application. Bosonic SPT phases are often constructed by using 
mathematical machinery such as group cohomology/cobordism and their physical properties are typically explained via 
gauge/gravity anomalies [HI |T71I^H HZ] . In this paper, we shall restrict our attentions to rather simple realizations 
with Z2 symmetry and study their “properties” in a systematic manner. (Generalization to Zjv symmetry, or arbitrary 
abelian symmetry, is also possible). The proposed models have a short-range entangled unique gapped ground state 
on a closed manifold which is protected by higher-form global symmetry. In our construction, charged excitations of 
different dimensionality may coexist and form a short-range entangled state which is protected by symmetry operators 
of different dimensionality. For instance, to provide some insight, we mention the existence of a non-trivial (5 -I- 1)- 
dimensional model protected by 0-form, 1-form and 2-form Z2 (8)Z2 ® Z2 symmetry where charged particles, loops and 
membranes form a short-range entangled vacuum. Kapustin and Thorngren provided a framework for constructing 
SPT phases protected by 0-form and 1-form symmetries by using 2-group. A complementary feature of our work 
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is that the model admits arbitrary g-form symmetries and its construction is simple. Also, some of our models do 
not seem to ht into the framework of the 2 -group construction despite the fact that these models possess 0 -form and 
1 -form symmetries only. 

Let us clarify what we mean by non-trivial SPT phases. We will say that a wavefunction is non-trivial if its 
preparation from a product state requires a large-depth local quantum circuit. By a local quantum circuit, we mean a 
unitary operation that can be implemented by a local Hamiltonian which may be time-dependent. By a large-depth, 
we mean that it takes a system-size dependent time to implement the required unitary operation. Ground states 
of (intrinsic) topologically ordered Hamiltonians are known to be non-trivial in this sense. Non-triviality of SPT 
wavefunctions can be dehned by restricting considerations to local quantum circuits which commute with imposed 
symmetries. Namely, each local component of the quantum circuit needs to commute with symmetry operators. 
To prove the non-triviality of the wavefunctions in the presence of symmetry, we employ three different arguments. 
First we consider the models with open boundaries and demonstrate that they have protected boundary modes 
which trivial models cannot possess. Namely, boundary mode may be gapless, spontaneous symmetry-breaking phase 
and/or topological phase depending on dimensionality of symmetry operators and spatial dimensions. While this 
type of arguments has been extensively used as simple diagnostic of non-triviality of SPT wavefunctions on the 
bulk, its connection to the circuit complexity of wavefunctions is rather indirect. Second we minimally couple the 
SPT Hamiltonians to generalized (higher-form) gauge fields. Namely, for g-form symmetry, the system is coupled 
to (g -I- l)-form gauge fields. We show that gauged versions of trivial and non-trivial SPT Hamiltonians belong to 
different topological phases, which implies that the original Hamiltonians (and wavefunctions) belong to different 
quantum phases in the presence of symmetries. This type of argument was originally proposed by Levin and Gu in 
the context of 0-form SPT phases, and can be converted into a fairly rigorous argument by formulating the gauging 
as a duality map. (See [5B] for instance). Third we shall construct a gapped domain wall in a (d -I- l)-dimensional 
topological phase by gauging the d-dimensional SPT wavefunction in (d -I- 1) space dimensions by following the idea 
in [29| . We then show that a gapped domain wall transposes excitations in a non-trivial manner which is possible 
only if the underlying SPT wavefunction is non-trivial. This type of argument is perhaps non-standard in condensed 
matter and high energy community, but is indeed natural from quantum information theoretical viewpoint. These 
three characterizations show that proposed models of higher-form SPT phases are indeed non-trivial in a sense of 
quantum circuit complexity. 

There are two key ingredients in our construction. First all the models are constructed on d-dimensional simplicial 
lattices which are d -I- 1-colorable, meaning that one can assign d -I- 1 distinct color labels to vertices of the lattice 
such that neighboring vertices have different colors. Golorable graphs have a number of useful properties in discussing 
lattice realizations of topological theories pTiUHH] . For instance, we shall see that gauged models can be defined on the 
same graph without modifying the lattice structure. The simplest example of colorable graphs is a two-dimensional 
triangular lattice where three color labels can be assigned to vertices. Second, to construct non-trivial models, we 
shall use unitary phase gates, known as multi-qubit control-Z gates in quantum information community, which are 
closely related to a certain non-trivial m-cocycle function for G = (Z 2 )®"* [221 IS]- Multi-qubit control-Z gates are 
of particular importance in the context of fault-tolerant quantum computation since they are outside of the Glifford 
group when involving three or more qubits |35j . 

The overall construction is summarized as follows. Let ai, 02 ,..., a^+i be d -I- 1 distinct color labels. We think of 
splitting color labels into several groups. We then place qubits on g-simplexes according to the splitting. For instance, 
for d = 5, one may have the following splitting: 


ai|a2a3|o4a5a6 


( 2 ) 


and qubits are placed on 0-simplexes of color oi, 1-simplexes of color 0203 and 2-simplexes of color 040505 . We then 
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apply multi-qubit control-^ gates on each d-simplex to obtain the non-trivial Hamiltonian from a trivial Hamiltonian. 
In the above example, we will obtain a five-dimensional model with (0,1, 2)-form Z 2 ® Z 2 ® Z 2 symmetry by applying 
three-qubit control-^ gates. The system supports point-like, loop-like and membrane-like charged excitations. While 
our studies are limited to models with Z 2 symmetries, the constructions can be generalized to systems with Z^v 
symmetries by using a certain generalization of control-^ gates as briefly explained in Section [V] Our construction 
can be viewed as a special realization of hyper-graph states recently proposed in |36j . and thus, applications to 
measurement-based quantum computations may be an interesting future problem. 

Classification of topological phases of matter is a problem of fundamental and practical importance, bridging 
condensed matter physics and quantum information science. At a formal level, lattice models of topological phases 
of matter can be probably classified by using the framework of higher-category theory. This, however, does not 
mean that classification of topological order is completed since category theoretical approaches provide only a set of 
consistency equations, such as pentagon and hexagon equations. Solving consistency equations are rather difficult 
both analytically and computationally, and thus finding a non-trivial solution to consistency equations seems to be 
the real challenge. Our lattice models, before and after coupling to gauge fields, presumably satisfy these consistency 
equations of category theoretical approaches, and are beyond known theories of topological order, such as the Walker- 
Wang model m and the Dijkgraaf-Witten topological gauge theories [T]. Thus, our model may serve as a stepping 
stone to further looking for exotic topological phases of matter, which may be of importance for quantum information 
processing purposes. 

Indeed, proposed models of SPT phases with higher-form symmetry have interesting quantum coding applications. 
In [55J [2^ , it has been pointed out that the classification of bosonic SPT phases with 0-form symmetry and classifica¬ 
tion of fault-tolerant logical gates in topological quantum codes are closely related. The construction of bosonic SPT 
phases with higher-form symmetry enables us to find (somewhat surprising) fault-tolerantly implementable logical 
gates. For instance, consider a four-dimensional system consisting of two copies of the (l,3)-toric code and one copy 
of the (2, 2)-toric code. Here, the (a, 6)-toric code refers to the (a + 6)-dimensional toric code with a-dimensional 
Pauli-Z logical operators and ^-dimensional Pauli-A logical operators. We choose these three copies of the toric code 
to be decoupled from each other. As we shall see, there exists a non-trivial three-dimensional SPT model which 
is protected by 0-form Z 2 ® Z 2 symmetry and I-form Z 2 symmetry. The presence of such an SPT phase implies 
that one can implement a three-qubit control-Z logical gate, belonging to the third-level of the Clifford hierarchy, 
among three copies of the four-dimensional toric code fault-tolerantly by a finite-depth local quantum circuit. This is 
rather surprising given the fact that the (1,3)-toric code and the (2, 2)-toric code possess logical operators of different 
dimensionality and thus belong to different topological phases. Our approach may give a hint on how to implement 
multi-qubit unitary logical gates on multiple quantum error-correcting codes of different code generators. 

The paper is organized as follows. Section [H] of the paper is devoted to a brief review of bosonic SPT phases with 
0-form symmetry. While this is a throughly explored subject, we try to provide a concise, yet precise treatment 
of various properties of 0-form SPT phases, such as the boundary mode, outcome of coupling to gauge fields and 
corresponding gapped boundaries, by using concepts from quantum information theory. Although we do not provide 
any new results in this section, our treatment and observations may have some novelty. In section [Till we will present 
lattice realizations of a bosonic SPT phase with generalized global symmetry and discuss their physical properties. 
We demonstrate that boundary of these models exhibits spontaneous breaking of symmetry, gapless critical modes 
or/and topological phases. In section IV we comment on quantum coding implications of our results. 



4 


II. TOPOLOGICAL PHASES WITH O-FORM SYMMETRY 

We begin by presenting a simple recipe of how to construct certain families of bosonic SPT wavefunctions with 
0-form symmetry by using multi-qubit control-Z operators. Namely, we construct wavefunctions for d-dimensional 
SPT phases with symmetry. We then study their boundary modes by finding dressed boundary operators 

which commute with the bulk Hamiltonian while satisfying certain commutation relations of Pauli operators. In a 
quantum coding language, boundary modes and dressed boundary operators can be viewed as a codeword space and 
logical operators of the bulk Hamiltonian respectively. This allows us to study symmetry constraints on boundary 
terms and discuss the boundary mode protected by symmetry. We also discuss the procedure of coupling SPT phases 
to gauge fields and corresponding gapped domain walls. 


A. Multi-qubit control-.^ gates 

We begin by presenting the definition of a multi-qubit control-^ gate [29l [34] (see Fig. [^a) for its quantum circuit 
representation). The control-^ gate, denoted by CZ, is a two-qubit phase gate acting in the computational basis as 

CZ\x,y) = {-iyy\x,y) x,y = 0,l (3) 

which adds —1 phase if both the first and second qubits are in |1) state. One can generalize the control-^ gate to a 
system of multiple qubits. The n-qubit control-Z gate, denoted by C^^~^Z, acts as follows 

C®”-'Z|a;i,...,x„) = (-l)"-""|cri,...,a;„) x, =0,1. (4) 



FIG. 1: (a) A quantum circuit representation of the three-qubit control-^ gate C^^Z. (b) Conjugation by the three-qubit 
control-.^ gate. Note that = C^'^Z. 

It is convenient to summarize how Pauli X operators transform under conjugation by multi-qubit control-Z oper¬ 
ators. For a two-qubit control-Z gate, one has 


CZ{Xi)CZ = {Xi)Z2 CZ{X2)CZ = Zi{X2). 


( 5 ) 
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For a multi-qubit control-Z gate, one has 


= (Xi)C®”"^Z2...„. (6) 

Here acts on the jth qubits {2 < j < n). So, conjugation by adds “decoration” of on 

Pauli X operators (Fig. [^b)). 

Multi-qubit control-Z gates have particularly useful applications in quantum coding theory since belongs 

to the n-th level of the so-called Clifford hierarchy [SniEllMlIMSS] which is an important concept in classifying fault- 
tolerantly implementable logical gates in topological stabilizer codes. Readers who are familiar with topological gauge 
theories may recognize the similarity between operators and a non-trivial n-cocycle function for G = Z 2 ®"': 

= (—l)®i'••■9""^ where = {g[^\ g 2 \ ■ ■ ■, gn^) and =0,1 [H SHI EH- For the connection 
between group cohomology and the Clifford hierarchy, see |28j . 


B. One-dimensional model with Z 2 ® ’^2 symmetry 

In this subsection, we study the one-dimensional SPT phase with Z 2 C) Z 2 symmetry [ 3111 . Consider a one¬ 
dimensional chain of 2n qubits with periodic boundary conditions. We assign color labels a, b to vertices in a bipartite 
manner such that odd (even) sites have color a (b). The Hamiltonian is given by 

2n 

= - E Oj, O, = X,.,Z,X,+,. (7) 

Since = 0, the ground state \tjj) satisfies Oj\ijj) = \ip) for all j. The Hamiltonian and the ground state have 

Z 2 0 Z 2 symmetry corresponding to two symmetry operators: 

n n 

(8) 

i=i 1=1 

where act on vertices of color a,b respectively. We can see that the Hamiltonian respects the symmetry: 

= 0, as well as the ground state: = IV')- To verify this, observe that 

= n”=i 02,-1, 5 w = nE02,- 

One can see that the ground state |V') is short-range entangled by considering the following finite-depth quantum 
circuit: 


eeE 


( 9 ) 


where CZe acts on two qubits on the edge e and E represents the set of all the edges. The superscript in f/O-O) 
indicates that we construct a model with two copies of 0-form Z 2 symmetry. One has 

= C7(o.o)^jj[/(o.o)t^ i7o = -^Xj, (10) 

j 

where |V’) = |-|-)®^” and |-|-) := ^(|0) -I- |1)). However, the quantum circuit is not symmetric since each 

local component, CZe, does not commute with symmetry operators 5'^“^ or (One important subtlety is that 

lj{o,o) commutes with and as a whole, but one cannot implement it through local quantum gates which 
commute with 5”^“^ and Indeed, one can verify that there is no finite-depth symmetric quantum circuit which 
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creates |'0) from a product state. See [35] for instance. In this sense, we say that |?/;) is a non-trivial SPT wavefunction 
in the presence of Z 2 0 Z 2 symmetry. 

We then study the boundary mode for the one-dimensional Z 2 0 Z 2 SPT phase by following an approach used by 
Levin and Gu m- Consider a one-dimensional chain of 2n spins with boundaries. Let bulk and boundary be sets 
of vertices in the bulk and on the boundary respectively as shown in Fig. We labeled qubits by oi, 61 ,..., On, fon, 
and boundary = {ai, 6 „}. One can write a generic form of the Hamiltonian as follows: 

^ ~ ^boundary T ^bulk ~ J ^ ^ Oj (H) 

jGbulk 

where Hboundary iuvolves terms localized near the boundaries and Oj are given by Eq. Q. We assume that the 
interaction strength J > 0 is sufficiently large so that one can restrict attentions to the low-energy subspace C: 

C = {\iIj) : Ojlip) = IV’) j G bulk} . (12) 

The low-energy subspace C is a 4-dimensional space. In the large J limit, one can assume that the boundary terms 
commute with the bulk Hamiltonian: 


{^boundary ^ ^hulk\ b- 


Later, we will discuss the cases where Hyoundary does not commute with H^uik- 


boundary bulk boundary 

ai bi an bn 



(13) 


FIG. 2: Bulk and boundary vertices on a one-dimensional chain with boundaries. 

We hope to find operators which characterize the boundary mode associated with the low-energy subspace C. 
Ordinary boundary Pauli operators, Xa ^, Za ^, Xy ^, on the boundary, do not commute with Htuik, and thus are not 
appropriate operators to describe physics of the boundary mode. One needs to find a complete set of Pauli operators 
which act non-trivially inside C, but preserve C. Let us consider the quantum circuit which is truncated at the 

boundary: = Usge = CZ^.b^CZy^a, ■ • • CZ^^.^a. 

the role of Pauli operators: 

One can verify that these Pauli operators commute with Oj for j G bulk. As such, Xa^, Za^ characterize the left 
boundary mode and Zy^ characterize the right boundary mode. In the language of quantum coding theory, one 
may view C as a codeword space, and Xai, Zai, Xb„, Zy^ as logical operators acting non-trivially inside C. In this 
picture, boundary modes are logical qubits encoded in the codeword space C. These boundary operators are “dressed” 
in a sense that they involve Pauli operators on the bulk. 

While dressed boundary operators commute with the bulk terms, they are not symmetric under 5”^“^ and 


CZa^b„- The following dressed boundary operators play 

J, 

Zb^ = = Zb^. 
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Indeed, dressed boundary operators are transformed under conjugation by 5'^“^ and as follows: 

^ ^ -Za, Zb„ ^ -Z6„ 

: Xai ^ai ^ai -^b„ 

From the above relations, one can deduce the action of symmetry operators inside C as follows: 

- Z,, 0 Z 6 „, ~ Z,, 0 Zb„. (16) 

The notation implies that two operators act in an identical manner inside the low-energy subspace C. We can 
then ask what kinds of boundary terms are allowed under this symmetry. Observe that a term on the left boundary 
needs to commute with Xa^ and Za ^, implying that there is no term that can be added on the left boundary. A similar 
argument holds for the right boundary. The only possible terms for Hboundary are and which are highly 
non-local. Therefore, the degeneracy on the edges cannot be lifted. One can consider boundary terms Hyoundary which 
do not commute with Hy^ik too. In such cases, perturbative analysis implies that non-trivial coupling between four 
degenerate ground states appear only in the 0 (n)-th order perturbative expansion which is exponentially suppressed. 
Therefore, one expects that the energy splitting among four low-energy states is exponentially small with respect 
to the system size n. The conclusion is that four degenerate boundary states are protected by symmetry and the 
degeneracy cannot be lifted by small perturbations which respect the imposed symmetry. 


Zb„ ^ ~Zbr, ■ 


(15) 


C. Two-dimensional model with Z 2 ® 2 , 2 ® 1-2 symmetry 

In this subsection, we shall study the two-dimensional SPT phase with Z 2 ( 8 )Z 2 ( 8 )Z 2 symmetry. Consider a triangular 
lattice as depicted in Fig. which is three-colorable in a sense that one can assign color labels a, b, c to vertices in 
such a way that neighboring vertices have different color labels. Qubits are placed on vertices of the lattice. Consider 
the trivial Hamiltonian Hq = — where V represents the set of all the vertices. We shall apply the following 

finite-depth quantum circuit, consisting of CCZ operators acting on each triple of qubits contained inside a triangle: 

17(°’0.0) ;= CCZ,,,.fc (17) 


where A represents the set of all the triangles. The resulting non-trivial Hamiltonian can be written as 

0„=A„ ll CZ, (18) 

vGV eGl-link(D) 

where 1 - link(u) represents the set of all the 1-links of a vertex v (Fig.[^. (A 1-link of v forms a 2-simplex by adding v. 
For a precise definition of 1-link, see [53]). Thus, interaction terms are Pauli X operators decorated by CZ operators. 
This Hamiltonian has Z 2 0 Z 2 0 Z 2 symmetry with respect to the following triple of Z 2 symmetry operators: 

(g) A„, (g) A„ (g A„ (19) 

'uGT<“) 

where represent the sets of all the vertices of color a, 6 , c respectively. That is, Pauli X operators act 

on qubits with distinct colors. Direct calculation shows that the Hamiltonian respects the Z 2 ® Z 2 0 Z 2 symmetry: 
[i7i, = [iJi, = 0. Also, the ground state l^/i) is symmetric since Ilweyf'’) ~ 

~ obtain this, recall that {CZ^ = I. As such, the Hamiltonian Hi possesses Z 2 0 Z 2 0 Z 2 



FIG. 3: Two-dimensional SPT phase with 1^2® 1^2® ^2 symmetry. Each of Z 2 symmetry operators acts on qubits of distinct 
colors. The quantum circuit for this system consists of CCZ acting on each triangle. The figure on the right shows an interaction 
term which is a Pauli X operator at a vertex v decorated by CZ operators acting on 1-links of v. 


symmetry. Under Z 2 0 Z 2 0 Z 2 symmetry, there are 128 = 2^ different SPT phases with 7 distinct generators which 
can be sorted into three types, called type-I, type-II and type-III. Upon gauging, type-I and type-II models are dual 
to abelian quantum double model with semionic statistics while type-III model is dual to the non-abelian ZI 4 quantum 
double model nsniiiiiiia. It has been shown that the aforementioned model corresponds to the SPT phase which 
contains the type-III cocycle function [55j. 

We then study the boundary mode for the two-dimensional SPT phases with Z 2 0 Z 2 0 Z 2 symmetry. Consider 
the triangular lattice with a boundary as shown in Fig. where boundary contains vertices of color a or & on the 
boundary. We write the Hamiltonian as H = H boundary + where consists of all the terms with v S bulk 

and [Hi,oundary^ Hy^ik] = 0. We construct dressed boundary operators which commute with the bulk terms and have 
proper commutation relations of Pauli operators. Let U be the quantum circuit which is truncated at the boundary 
such that U is a product of CCZ acting on all the triangles which are dully contained on the lattice. By applying 
the quantum circuit U to Pauli operators on edge qubits, one can construct dressed boundary operators. Pauli X 
operators are decorated with control-Z operators involving bulk qubits as shown in Fig. |^b) while Pauli Z operators 
remain unchanged. We shall denote these dressed boundary operators by Xa ^, Xi,j , Za^, Zy-. 

Let us study how dressed boundary operators transform under symmetry operators, S^°'\ : 


Sia) 

sw 

sC) 


X, 


X 


Xy. 

a,- 5 X^bj 


Xa, ^ X 

^aj Zf, Xa.Zi, 


Xb,, 

Xb„ 


Zaj —Zaj, 

Zaj Za^ , 

Zaj Z^bj Z j 


From these relations, one hnds 


Zbj Zbj 
Zb^ -t -Zbj 
Zax Za., 




( 20 ) 


( 21 ) 


eGE 


where E represents the set of all the edges on the boundary and the products for S^°‘\ run over all the vertices 
of color a, b on the boundary. Let us write down possible boundary terms which respect the symmetry: 


Zaj Zaj 


Zb, Zh 


Zaj + Zhj_xXajZbjj Xbj+ZajZyX, 


“j + i ■ 


( 22 ) 
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FIG. 4: (a) Edge and bulk qubits, (b) Dressed boundary operators Xa and Xi,. (c) Dressed boundary operators after 
conjugation by a symmetry operator 

The first two terms are ferromagnetic interactions among qubits of color a or 6 respectively. The third and the last 
terms lead to the following Hamiltonian at quantum criticality: 

2n 2n 

H = -J2 Zj-iX,Z,+, -Y.X,. (23) 

t=i i=i 

This Hamiltonian can be transformed into two decoupled copies of critical quantum Ising model by a duality trans¬ 
formation Therefore, one can conclude that the boundary mode may support Z 2 ferromagnets (i.e. spontaneous 
breaking of Z 2 0 Z 2 ), or, Z 2 critical models with gapless modes. It is interesting to observe that the above critical 
Hamiltonian can be written as a sum of one-dimensional SPT Hamiltonians with Z 2 ® Z 2 symmetry: H = Hq + Hi 
where Hq = - Xj and Hi = - J2^jZi Zj-iXjZ^+i and the symmetry operator is identical to the quantum 

circuit for a non-trivial Z 2 0 Z 2 SPT phase: = YieeE Namely, one has H = Hq + 


D. Coupling to gauge fields 

In this subsection and the next, we study the procedure of “gauging” m- Physically, gauging is a process of 
minimally coupling a system with global symmetry G to gauge fields with gauge symmetry G. Formally, gauging 
can be viewed as an isometric bijective map (ie. a duality map) from wavefunctions with global symmetry to 


^ Consider the following transformation: Zj-iXjZjj^i —> Zj — iZjj^i and Xj —> Xj. Strictly speaking, this duality transformation is well 
defined only for an infinite system. 
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wavefunctions with gauge symmetry m- Detailed discussions on physical properties of gauged models are beyond 
the scope of the present paper. Instead, we will present a generic framework of gauging non-trivial models of bosonic 
SPT phases. We shall see that the use of colorable graphs is particularly convenient in constructing gauged models. 

In this subsection, we briefly review the procedure of Z 2 gauging. Consider a square lattice on a torus where qubits 
live on vertices, and consider the following trivial system with Z 2 symmetry: 

H = -Y,^v (24) 

V 

where v represents vertices. Clearly, this Hamiltonian and its ground state is symmetric under Z 2 global transformation 
S = Xy. By applying the Z 2 gauging map, one can transform this trivial system with global Z 2 symmetry into a 
two-dimensional system with Z 2 gauge symmetry, namely the two-dimensional toric code, as shown below. 


a b (b) a+b 

O-O -O- 



FIG. 5: Qubits and interaction terms before and after the Z 2 gauging. 

To begin, consider a system of two qubits \a, b) where qubits live on endpoints (vertices) of an edge and a, 6 = 0,1 S 
Z 2 (Fig.§. Imagine a transformation (C^)®^ —whose output is a one-qubit state |a -I- &) where the summation 
is modulo 2 and the output qubit lives on the edge (Fig. [^. Next, consider a system of qubits supported on vertices 
of a square lattice and denote the entire Hilbert space by "Ho- Consider a computational basis state, | 5 i, 52 , • ■ • )5ra) 
where Qj = 0,1, in "Ho- We think of applying the above transformation to every edge of the system and obtaining 
an output wavefunction living on edges of the square lattice. This defines a map F from computational basis states 
in 'Ho to computational basis states in "Hi where T-Li represents the Hilbert space for a system with qubits living on 
edges of a square lattice. We have 


F(|5 i; • ■ ■) 5 n)) — |hi,..., hn'), Sj,hj — 0,1 (25) 

where n, n' are the numbers of vertices and edges respectively. An important property of F is that the output 
wavefunction satisfies the gauge constraints. Namely, for an arbitrary plaquette operator Bp — n y^pZy, one has 
I ^1 5 ■ ■ ■ 5 ) - I hi , ... , hyt ) . 

One can extend the gauging map F to wavefunctions which are not computational basis states. For this purpose. 
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we will consider some subspaces of "Ho and "Hi. Define the symmetric subspace of Hq as 


= m e Ho : S\^) = IV')} (26) 

where S = 0^, Xy. Define the gauge symmetric subspace of Hi as 

Hf™ = {|^) G Hi : B{^m = \i^) Vy} (27) 

where 7 represents an arbitrary closed loop on a square lattice and - 6 ( 7 ) = rije 7 ^3 i® ^ Wilson loop operator. Here 
we consider not only contractible loops, but also arbitrary closed loops which may have non-trivial winding. We then 
define the gauging map as follows 

r(IV')) = 4 E Cg,_gSi\9i,---,9n)) (28) 

where |V') = J2g,,...,g^ C'si,..., 3 „ Isi, ■■■,9n) & Hq*"". Then one has 

dim Ho^"* = dim Hf . (29) 


Even more, the gauging map T is bijective and isometric {i.e. is a duality map). By an isometric map, we mean that 
the inner product of any pair of wavefunctions is preserved. 

Let us apply the gauging map to a trivial wavefunction \tp) of Eq. (24). Let |i/>) be the output wavefunction. We 
shall see that Ii/j) is a ground state of the two-dimensional toric code: H = — ~ J2p ^p- Indeed, Z 2 gauge 

constraints account for the plaquette terms Bp = UeepZe in the toric code. As for the star terms Ay = UyeeXe, 
observe that Xylijj) = \tp) inside Hg*^™. Flipping a spin at the vertex v is equivalent to flipping four spins on edges 
that are connected to the vertex v in the gauge theory. Thus, Xy operator in Hg*^™ is equivalent to the vertex term 
Ay in . As such, the output wavefunction must satisfy Ay\ip) = \ip) for all u, and thus is a ground state of the 
toric code. The above procedure can be extended to a d-dimensional system with on-site symmetry group G where G 
is an arbitrary finite group, and the gauging map outputs the d-dimensional quantum double model with G [Hill]- 


E. Gauged model and gapped domain wall 

In this subsection, we apply the gauging map defined in the previous subsection to SPT wavefunctions. Consider 
the two-dimensional SPT phase with Z 2 0 Z 2 0 Z 2 symmetry supported on a three-colorable lattice where qubits 
are placed on a,b,c vertices respectively. Let us denote the entire Hilbert space by Hg. The gauging map can be 
expressed as follows 


P = Pa 0 Pb (g) Pc (30) 

where r(j,r;,,rc are Z 2 gauging maps acting on qubits living on a, 6 , c respectively. Recall that edges of a colorable 
graph can be labeled by pairs of color indices, ab, be, ca, by looking at colors of vertices that are connected by edges. 
Imagine that we place qubits on ab, be, ca edges, instead of a, b, c vertices, and denote the entire Hilbert space by 
Hi (Fig. [^a)). Observe that input wavefunctions of the gauging map Tq live on vertices of color a while its output 
wavefunctions live on edges of color be since the middle point of two neighboring vertices of color a is an edge of 
color be (Fig. |^b)). Thus, the gauging map P is a map from computational basis states in Ho to those in Hi. A 
key observation is that, due to the colorability of the graph, one does not need to modify the lattice structure in 
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defining the Hilbert space "Hi for the output wavefunctions. One can define Z 2 0 Z 2 0 Z 2 symmetric subspace 
and Z 2 0 Z 2 0 Z 2 gauge symmetric subspace as before. Then the gauging map T = Tq 0 T;, 0 Tc is a duality 

1 , n/sym 1 n/sym 

map between Tiq and . 


(a) (b) 



FIG. 6: (a) Qubits on edges of color be. (b) Gauging maps Fa. (c) Interaction terms. 


Let us apply the gauging map F to the two-dimensional SPT phase. Let IV’) be a ground state of the two-dimensional 
SPT Hamiltonian with Z 20 Z 20 Z 2 symmetry and |'0) = r(|'0)) be the output state. Let be a plaquette surrounded 
by vertices of color a, as shown in Fig. |^a), and be the set of such plaquettes. Then one has i?p(a)|(/)) = \(j)) 
where Bp{a) is a tensor product of Pauli Z operators surrounding the plaquette . Similar argument holds for Bp{b) 
and Bp(o). The original symmetric wavefunction |i/i) satisfies X.^^a)\ip) = X.^(b}\4’) = X.^(a)\'4’) = I'*/')- After gauging, 
one can write down the corresponding vertex terms: 

A^(a) l-tp) = A^(b) l-tp) = A„(c) 1'!^) = I'lp). (31) 

which resemble ordinary vertex terms for the toric code, but are additionally decorated by CZ operators as shown in 
Fig. I^c). As such, the gauged model can be written as 

H = - ^ ^P~ A„(a) - ^ ^ A„(c). (32) 

pgp{a)^p{t)_p(c) i,(<i)gy(a) t,(li)gy(i>) i,(c)gy{c) 

The Hamiltonian can be viewed as three copies of the toric code which are intricately coupled with each other via CZ 
phase operators. 

An interesting application of the gauging map is that one can construct a gapped domain wall in a d-dimensional 
topologically ordered system by using bosonic {d — l)-dimensional SPT phases [551 HH]- To be specific, consider a 
two-dimensional colorable graph A with color labels a, 6 , c where qubits live only on vertices of color a, b. Consider 
a one-dimensional line dA in the graph which consist only of vertices of color a, b as shown in Fig. [^a), which splits 
the entire system into the upper and lower parts. Imagine that a one-dimensional SPT wavefunction with Z 2 0 Z 2 
symmetry lives on i9A while all the other qubits are in the trivial product state of |-|-). 

We couple the entire system to Z 2 0 Z 2 gauge fields living in two dimensions. In a gauge theory, qubits live on 
6 c-edges and ca-edges where the gauging map F acts on vertices of color a, b, and outputs quantum states on 6 c-edges 
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FIG. 7: (a) A one-dimensional SPT wavefunction in two dimensions and the corresponding gapped domain wall, (b) The 
gauging map. (c) Terms near the domain wall. 

and ca-edges as shown in Fig.j^b). The Hamiltonian for the gauged model can be written as 

H = Hup + Hdown + Hqa (33) 

where Hup, Hdown represent terms on A\dA while Hqa represents terms connecting two Hamiltonians Hup and Hdown- 
Note that Hup, Hdown are identical to those of two copies of the toric code while Hqa is different from those of the 
ordinary toric code and can be viewed as a domain wall connecting Hup and Hdown- Terms in Hqa are explicitly 
written in Fig. [^c). Let Ca, Cb be electric charges associated with violations of vertex-like Pauli-X terms. Let iJia, mb 
be magnetic fluxes associated with violations of plaquette-like Pauli-Z terms. Then labels of anyonic excitations get 
transformed upon crossing the domain wall as follows: 

ma -t maCb, mb -t mbCa, -)► e^, eb -t Cb- (34) 

The fact that the domain wall transposes labels of anyonic excitations implies that the one-dimensional SPT 
wavefunction is non-trivial and cannot be created from a trivial state by a symmetric finite-depth quantum circuit [29| . 
The key observation is that, while the domain wall is localized along a one-dimensional region, it cannot be created 
by a local unitary transformation acting on qubits in the neighborhood of the domain wall. Suppose there exists a 
local unitary U which creates the domain wall by acting only on qubits in the neighborhood of the domain wall. Let 
£ be a string operator corresponding to the propagation of a magnetic flux in the absence of the domain wall. Then 
differs from i only at the intersection with the domain wall. This implies that a magnetic flux remains to be a 
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magnetic flux upon crossing the domain wall, leading to a contradiction. Thus, to create the domain wall, one needs 
to apply a local unitary transformation on all the qubits on one side of the system. This argument enables us to 
show that the underlying one-dimensional SPT wavefunction is non-trivial. Suppose that there exists a symmetric 
local unitary transformation U which creates the underlying SPT wavefunction. After gauging this symmetric unitary 
operator, one obtains local unitary transformation U defined in the gauge theory which is localized along the domain 
wall, leading to a contradiction. In the above argument, braiding statistics of anyons plays the role of topological 
invariants in proving the non-triviality of SPT wavefunctions. 

Finally, we briefly comment on a gapped domain wall which can be constructed by gauging the two-dimensional SPT 
wavefunction in three dimensions. The domain wall connects three copies of the three-dimensional toric code where 
electric charges ei, 62,63 are point-like while magnetic fluxes mi, m 2 , m 3 are loop-like. Upon crossing the domain wall, 
electric charges remain unchanged: 


ei ei 62 62 63 63 (35) 

while magnetic fluxes transform into composites of magnetic-fluxes and loop-like superpositions of electric charges: 

mi —>■ miS23 m 2 —>■ m2Si3 W 3 —>■ m3Si2- (36) 

Here Sij (i ^ j) is a one-dimensional excitation which is a superposition of electric charges et and Cj. It has been found 
that Sij can be characterized by a non-trivial wavefunction of a one-dimensional Z 2 (E) Z 2 bosonic SPT phase [29] . 
Namely, if the emerging wavefunction is written as a superposition of excited eigenstates, its expression is identical to 
a fixed-point wavefunction of a one-dimensional non-trivial SPT phase. It is important to note that these fluctuating 
charges sij are loop-like objects which are unbreakable in a sense that their creation requires membrane-like operators. 
In addition, m^ and sjk exhibit non-trivial three-loop braiding statistics 


F. Higher-dimensional generalization 


Finally, we briefly present construction of d-dimensional SPT phases with Z 2 ®'^~''^ symmetry. Consider a d- 
dimensional simplicial lattice which is {d+ l)-colorable with color labels oi,..., ad+i and place qubits on each vertex. 
Let be a trivial Hamiltonian. The non-trivial Hamiltonian Hi can be constructed from the 

following finite-depth quantum circuit: 


H 


(0,0,...) ^ jj(0.0,...)^(0,0,...)^(0,0....)t ^(0,0....) _ 


n 


C®‘^Zi 




(37) 


(ii,Z2,...,id+l)6A 


where A represents the set of all the d-simplexes. The interaction terms of are 


^( 0 . 0 ....) 


Oy=Xy [] 


(38) 


where {d — l)-link(u) represents the set of all the (d — l)-links of the vertex v. Namely, a (d — l)-link of u is a 
(d — l)-simplex which forms a d-complex by adding the vertex v. There are d -I- 1 copies of Z 2 symmetry operators 


^ In a three-loop braiding, two loops are braided while the third loop pierces through two loops m- 
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associated with color labels aj: 

n n ( 39 ) 


where represents the set of all the vertices of color aj. 

Boundary mode in higher-dimensional SPT phases can be studied in a similar manner. Namely, the (d — 1)- 
dimensional boundary can support a quantum critical Hamiltonian of the form H = where 

and are trivial and non-trivial d— 1-dimensional SPT Hamiltonians with symmetry respectively. The 

gauged model can be defined on the same colorable lattice where qubits are placed on centers of {d — l)-simplexes 
instead of vertices. The model looks like d-l-1 copies of the toric code whose vertex terms have decorations of 
operators, mixing d -|- 1 copies in an intricate way. If one gauges the d-dimensional SPT phase in d -|- 1 dimensions, 
one obtains a gapped domain wall in d-l-1 copies of the toric code where codimension-1 magnetic flux get transformed 
into a composite of codimension-1 magnetic flux and codimension-1 fluctuating superpositions of electric charge upon 
crossing the domain wall |29j . 


III. TOPOLOGICAL PHASES WITH GENERALIZED GLOBAL SYMMETRIES 

In this section, we present examples of bosonic SPT phases with higher-form global symmetries. The key distinction 
between models with 0-form symmetries and higher-form symmetries is that qubits are placed on g-simplexes for 
models with g-form symmetries. We use the multi-qubit control-Z gate to construct non-trivial SPT wavefunctions 
with g-form symmetries. We also study the boundary mode, gauged models and corresponding gapped domain walls. 


A. Three-dimensional model with 1-form symmetries 

In this subsection, we present a model of a three-dimensional SPT phase with 1-form Z 2 0Z2 symmetry. Consider a 
three-dimensional simplicial lattice A which is four-colorable with color labels a, b, c, d. Edges of the graph are labeled 
by pairs of colors, such as ab, ac, ad ,.... We place qubits on a6-edges and cd-edges of A while no qubits are placed on 
edges of other colors or vertices as shown in Fig. |^a). 

We specify global symmetry operators which have geometries of two-dimensional closed manifold. Consider a closed 
2-manifold Ad which intersects with edges of A. One may view AI as a two-dimensional simplicial sublattice of A. 
We place Ai such that it does not intersect with vertices. (In other words. Ad is a discretization of a plane on a dual 
lattice). Let 5(Ad) be a sheet of Pauli X operators acting on edges which are intersected by Ad. The simplest example 
of such a global symmetry operator can be constructed by considering a small 2-sphere as shown in Fig. [^b). There 
are two different types of symmetry operators, acting on ad-edges and cd-edges respectively, and they are separable. 
Namely, one can choose Ad such that it does not intersect with cd-edges. One can also glue 2-spheres, which are 
centered at vertices of color a and b, together to construct a sheet which intersects only with ad-edges. So the system 
has 1-form Z 2 0 Z 2 symmetry. 

We then construct the model Hamiltonian First the trivial Hamiltonian is given by Hq = where E 

represents the set of all the edges. The non-trivial Hamiltonian is obtained by applying CZ operators on every pair 


^ Our model seems very similar to the model proposed in US], which was treated as a bosonic 0-form SPT Hamiltonian protected by 
time-reversal symmetry. 
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(a) d (b) 




(c) (d) 



d 



FIG. 8: Three-dimensional SPT phase with 1-form Z 2 (8) Z 2 symmetry, (a) Qubits on a colorable graph, (b) A symmetry 
operator associated with a closed sphere, (c) A finite-depth quantum circuit, (d) Interaction terms. A Pauli-X is decorated 
by Pauli-^ operators acting on 1-links. 


of qubits inside each 3-simplex. Let us denote such a unitary operator by 

Yl CZ,,, (40) 

(ij)6A 

where A represents the set of all the 3-simplexes and CZ^ ^ acts on two qubits inside a 3-simplex (Fig. [^c)). The 
non-trivial Hamiltonian Hi = is given by 

Hi=-Y,Oe, Oe=Xe [] Z,, (41) 

e^E e'Gl-link(e) 

where l-link(e) represents the set of 1-links of an edge e (Fig. |^d)). Since is a finite-depth quantum circuit, 

Hq and Hi belong to the same topological phase in the absence of symmetry. Yet, each local component in 
does not commute with global symmetry operators. We claim that ground states of Hq and Hi cannot be connected 
by a symmetric local unitary. 

One can see that interaction terms Oe commute with 1-form symmetry operator S{M.) for any closed 2-manifold 
M. Let us verify that the ground state \tpi) of Hi is symmetric: S{M)\'4’i) = IV'i) for all M. We prove this for the 
cases where A4 is a contractible sphere. Let Aiy be a small 2-sphere surrounding the vertex v. With some speculation, 
one can confirm that 


\{Oe = S{M,) (42) 

v^e 

since all the Pauli Z operators cancel with each other. Since Oe|'0i) = IV'i)) one has S{M.v)\'4’i) = IV'i)- Any 
contractible sphere M. can be constructed by attaching M.v for various u, and thus 5'(AI)|V’i) = IV'i)- 

Under 1-form global symmetry, excitations are loop-like objects. Namely, to create excitations by operators which 
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commute with global symmetry operators, one needs to consider closed strings of Pauli Z operators. Let 7 “^ be 
a closed loop consisting of ab-edges and let ^( 7 “^) be a string of Pauli Z operators acting on 7 “^. This operator 
creates string-like excitations, violating Oe along 7 “^ while commuting with symmetry operators. A similar operator 
^(ycd) defined for a closed-loop of cd-edges. The 1-form symmetry, imposed by S{M), can be viewed as 

a conservation law for loop-like excitations where the number of the cuts of loop-like excitations made by A4 must 
be even. As such, there are two copies of Z 2 conservation law on loop-like excitations supported on ab-edges and 
cd-edges. 

We then study the boundary mode of the aforementioned three-dimensional model. We choose the two-dimensional 
boundary which consists only of vertices of color a, 6 , c such that the boundary can be viewed as a three-colorable 
graph with qubits living on ad-edges (Fig. [^a)). As before, we are interested in the Hamiltonian with the boundary 
term of the form: H = Htuik + Hboundary where Hboundary] = 0, and the low-energy subspace is denoted by 

C. Dressed boundary operators, denoted by Xg(ab), Z^(ab) for ad-edge can be found by applying the truncated 

quantum circuit U to Pauli operators associated with qubits on the boundary. We find that X^(ab) is decorated with 
a Pauli Z operator on the bulk while Z^(ab) = Z^(ab) remains unchanged (Fig. [^b)). 



(b) 


(c) 


c 



a 


FIG. 9: (a) Qubits on the boundary, (b) Dressed boundary operator. A vertex of color d is contained inside the bulk while 
vertices of color a, b, c are on the boundary, (c) Symmetry operators on the boundary. X and Z indicate that these operators 
are expressed in terms of dressed boundary operators. 


Let us study how symmetry operators 5'(rab) and S'(rcd) act on dressed boundary operators. Let us take Fab to 
be a sphere whose center is at the vertex v of color a or 6 on the boundary. Then S'(rab) acts as a vertex-like term as 
depicted in Fig. [^c). As for S'(rcd), consider a sphere with the center being at a vertex of color c on the boundary. 
Then is a plaquette-like term as shown in Fig. [^c). Thus the boundary mode supports topologically ordered 
states, namely the two-dimensional toric code while the bulk is trivial. The non-triviality of the boundary mode 
is a strong indication of the non-triviality of the three-dimensional SPT wavefunction. We shall present additional 
supports for the non-triviality below. 


B. Gauging 1-form symmetries 


We have seen that one can couple 0-form SPT phases to gauge fields where physical degrees of freedom lives on 
edges of a graph with gauge constraints acting on plaquettes. In the cases of 1-form SPT phases, one needs to modify 


gauge constraints as shown in Fig. 10 a) where physical degrees of freedom live on plaquettes and gauge constraints 


act on volumes. In this subsection, we illustrate the procedure of gauging 1-form symmetries. 
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Let us define the gauging map T for systems with 1-form Z 2 symmetries. To begin, consider a system of 12 
qubits which live on edges of a single cube as shown in Fig. [Io]^a). Consider a computational basis state of the form 
|a;i, 3 : 2 ,..., a;i 2 ) where Xj = 0,1. The output state of the gauging map T is a 6-qubit state where qubits live on 
plaqnettes of the cube. Namely, the spin values of the output state are given by Z 2 summation of spin values on edges 


surrounding the plaquettes as shown in Fig. 10 a). One can extend the above map to arbitrary three-dimensional 


lattices. Let Hi and H 2 be the Hilbert spaces for the systems where qubits live on edges and plaquettes of the cubic 
lattice respectively. Then F can be viewed as a map from computational basis states in 'Hi to those in H 2 . A key 
observation is that the output wavefunctions always satisfy generalized gauge constraints. Namely, if one sums up 
all the spin values on plaquettes surrounding a cube, then one obtains 0 modulo 2, implying Z 2 generalized gauge 
symmetry as in Fig. [lot^c). More generically, let A/” be a closed 2-manifold consisting of plaquettes of the lattice. 
We define the generalized gauge symmetry operator by T(A/’) which implements Pauli Z operators on plaquettes 
contained in JV. Then the output wavefunction |'0) always satisfies T{JV)\ip) = \tjj) for all Af. 

One can view F as a duality map by restricting our attentions to some proper subspaces of Hi and H 2 . Let 
Hi^™ C Hi be a subspace of wavefunctions which are symmetric under 1-form symmetry: 


Hf'" = {\iP) G Hi : S'(7W)|V’) = IV') VM}. 


(43) 


Here 1-form symmetry operators can be written as S{AA) for an arbitrary closed 2-manifold AA where S{AA) imple¬ 
ments Pauli X operators on edges intersected by M. (Fig. [Io|(b)). Similarly, we define the gauge symmetric subspace 
as follows: 


H^"’” = {|^) G H 2 : = 1^) VAA}. (44) 


Then one can see that F induces a duality map between H*^™ 
defined for the dual and direct lattices respectively. 




and H 2 ^”^. Note that closed 2-manifold Ad,A/” are 


\xi -\- X2 -\- 3 ^ 4 ) 




FIG. 10: (a) Gauging the 1-form Z 2 symmetry, (b) A 1-form symmetry operator, (c) A generalized gauge constraint. 

Let us apply the 1-form Z 2 gauging map to the trivial three-dimensional Hamiltonian: H = — Xe with a trivial 
ground state |'0) = |-|-)®”. The system respects 1-form Z 2 symmetry since the Hamiltonian consists only of Pauli 
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X operators. The original ground state |'0) satisfies Xelip) = \ijj) for all e. Upon gauging, the output wavefunction 
lip) satisfies Ae\'^) = |'0) where is a tensor product of Pauli X operators acting on plaquettes attached to an 
edge e. Thus the output wavefunction |?/;) is a ground state of the (2, l)-toric code. Recall that, if one gauges the 
trivial three-dimensional Hamiltonian as a model with 0-form symmetry, one obtains the (l,2)-toric code. In higher 
dimensions, the gauged models for the trivial Hamiltonian with 0-form and 1-form symmetry are the (1, d — l)-toric 
code and the (2, d — 2)-toric code respectively. 

Now, we shall apply the gauging map to the non-trivial three-dimensional Hamiltonian Hi with 1-form symmetry. 
One can write the gauging map as T = O Ted. Consider a plaquette consisting of a6-edges as shown in Fig. 11 'a). 
Let |a:i,..., X 2 n) be a computational basis state supported on ab-edges on the plaquette. The gauging map Tq;,, when 
acting on these qubits, outputs a single-qubit state \xi + X 2 + ■ ■ ■ + X 2 n) that lives at the center of the plaquette 
where summation is modulo 2. Notice that the center of a plaquette of cd-edges coincides with the center of a ab-edge 
as shown in Fig. [^a). Thus, for qubits supported on cd-edges, the output wavefunction is supported on qubits on 
ab-edges, and vice versa. Note that one does not need to modify the lattice structure in order to construct the gauged 
model. 



Let us find the gauged model of Hi. The original ground state wavefunction lipi) satisfies Oe\ipi) = \ipi)- The 
corresponding operators Of, are shown in Fig. 0b) which are identical to the original operators after exchanging 
Pauli X and Z operators. Namely, the gauged model can be written as 

Hi = UhHiUI (45) 

where Uh represents a transversal Hadamard operator. This implies that the gauged model Hi is not topologically 
ordered. Since Hq is topologically ordered while Hi is not, Hq and Hi belong to different quantum phases in the 
presence of 1-form symmetry. 

Next, we shall gauge the non-trivial three-dimensional Hamiltonian Hi in four dimensions and construct a gapped 
domain wall for the four-dimensional toric code. Consider a four-dimensional five-colorable graph with color labels 
a, b, c, d, e where qubits are placed on ab-edges and cd-edges. Choose a codimension-1 hypersurface (a three-dimensional 
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volume) dA and place the non-trivial 1-form SPT wavefunction with Z 2 0 Z 2 symmetry on dA while the rest of qubits 
are in the trivial state |-|-). We then gauge the entire system by coupling it to Z 2 0 Z 2 generalized gauge fields. To 
construct the output Hilbert space, we place qubits on centers of cde-plaquettes and a&e-plaquettes. The gauged 
model can be written as H = H^p + Hdown + H^aii where Hup and Hdown are identical to those of two copies of the 
four-dimensional (2,2) toric code while Hyjaii can be viewed as a gapped domain wall inserted along dA. 

Recall that the four-dimensional (2, 2) toric code has pairs of two-dimensional logical operators and possesses loop¬ 
like excitations. Let 61,62 and mi, m 2 be loop-like electric charges and loop-like magnetic fluxes in two copies of the 
four-dimensional toric code. Then one can verify that the domain wall transposes loop-like anyonic excitations as 
follows: 


ma -)■ maCb, mb -)■ mt,6a, 6a -)■ 6a, 6b 6b. (46) 

The non-triviality of the domain wall indicates that the underlying SPT wavefunction is non-trivial in the presence 
of 1 -form symmetry. 


C. Two-dimensional model with 0- and 1-form symmetry 

So far, we have considered SPT Hamiltonians where charged excitations are either point-like (0-form symmetry) 
or loop-like (1-form symmetry). It turns out one can construct a non-trivial SPT Hamiltonian protected by both 
0-form and 1-form symmetries where the system supports both point-like and loop-like charged excitations. In this 
subsection, we present a two-dimensional SPT Hamiltonian with 0 and 1-form Z 2 0 Z 2 symmetry. 

Consider a two-dimensional three-colorable graph with color labels a, b, c. We place qubits on ab-edges and c-vertices 
as in Fig. 12 'a). We start from the trivial Hamiltonian Hq = — ~ '"^here represent 

the sets of all the ab-edges and c-vertices respectively. We shall apply the following finite-depth quantum circuit: 


C/(04)= Y[ CZ,,j 


(47) 


(hi)6A 


where CZij acts on a pair of qubits contained inside a 2-simplex. The non-trivial SPT Hamiltonian is 


Hi=- Oe- ^ Oa. 

Interaction terms are Pauli X operators with some decorations of Pauli Z operators: 


(48) 


Oe=Xe Jl Zu 

•i;G0-link(e) 


Oy — Xy Z^(ab) 

eG l-link(i;) 


(49) 


where O-link(e), l-link(u) are sets of 0-links of e and 1-links of v respectively (see Fig. 12). 


The Hamiltonian has 0-form Z 2 symmetry and 1-form Z 2 symmetry acting on c-vertices and ab-edges respectively. 
The global 0-form symmetry operator is given by = JlaGyCc) where represents the set of all the vertices of 
color c. A 1-form symmetry operator is a string of Pauli X operators where A4 is a closed loop intersecting 

with ab-edges, S'^“*’^(A1) = IleGAi shown in Fig. 12'a). To see this, observe that S'^“*’^(A1) can be constructed 


by taking a product of Oe for edges crossed by M or contained inside M. Excitations associated with violations of 
Oe have geometries of closed loops due to the 1 -form symmetry. 

To verify the non-triviality of the SPT Hamiltonian, we study its boundary mode. Consider a one-dimensional 
boundary which consists only of vertices of color a,b as in Fig. ©a). See the figure for dressed boundary operators 
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(b) 



Z 


0-link 





FIG. 12: SPT phase with 0- and l-form symmetry, (a) A 1-form symmetry operator which is a closed string of Pauli X 
operators, (b) Interaction terms. Decorations are on 1-links and 0-links. 
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FIG. 13: (a) A boundary with a, 6-vertices, (b) A boundary with a, c-vertices. 


and the action of symmetry operators inside the low-energy subspace. Let us denote a6-edges on the boundary by ej. 
Allowed boundary terms are and their products, implying that the boundary mode supports a ferromagnetic 

order in the X basis. One can also consider a one-dimensional boundary which consists only of vertices of color a, c 
as in Fig. 17|b). Let us denote c-vertices on the boundary by Vj. Allowed boundary terms are Zy.Zy.^-^ which 


corresponds to a ferromagnet in the Z basis. Thus, in either choice of boundaries, the boundary mode would be a 
classical ferromagnet. Finally, one may consider a one-dimensional boundary where the aforementioned two types 
of boundaries coexist as depicted in Fig. See the figure for the action of symmetry operators in the low-energy 
subspace. The boundary mode, associated with a Xe^Xe^+i ferromagnet, is smoothly connect to the boundary 
mode, associated with a ZvXvj^i ferromagnet. Thus, properties of the boundary mode do not depend on choices of 
boundaries. 

Next, let us apply the gauging map. Since symmetry operators have different dimensionality, one needs to couple 
the system to gauge fields of different forms. Specifically, we shall consider the gauging map F = which 
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FIG. 14: A generic boundary consisting of a, b, c-vertices. The boundary mode consists of ferromagnets in A-basis and Z-basis 
which are smoothly connected. 

is a duality map between the symmetric subspace to the gauge symmetric subspace. Due to the colorability of the 
graph, the output of lives on c-vertices while the output of lives on a6-vertices. Let us gauge the trivial 
Hamiltonian Hq and non-trivial Hamiltonian Hi. The gauged model Hq consists of the toric code living on ab-edges 
and an Ising ferromagnet (the output wavefunction is the GHZ state) living on c-vertices. The gauged model Hi is 
identical to the original Hamiltonian Hi up to transversal Hadamard transformations. Gauged wavefunctions |^/’o) 
and \tjji) are not connected by a finite-depth quantum circuit, so Hq and Hi belong to different topological phases in 
the presence of symmetries. 

By gauging this non-trivial two-dimensional wavefunction in three dimensions, one can construct a non-trivial 
gapped domain wall. The system consists of the (1, 2)-toric code and the (2, l)-toric code with a domain wall. Anyonic 
excitations are given by eioop^epointi’m-ioop^'iTT-point- Upon crossing the domain wall, labels of anyonic excitations are 
transposed as follows: 


^Loop ^ ^loop ^point ^ ^point ^loop ^ ‘^loop^loop ‘^point ^ ^point^point • 


(50) 


D. Three-dimensional model with 0-, 0- and l-form symmetry 

In this subsection, we consider a three-dimensional SPT Hamiltonian with 0-, 0- and 1-form symmetries. 

Given a three-dimensional four-colorable graph with color labels a, 6, c, d, we place qubits on a-vertices, 6-vertices and 
cd-edges. Starting from a trivial Hamiltonian Hq = — X]uGy(“) ~ ~ non-trivial model 

is constructed by applying = n(i j fc)GA to the trivial Hamiltonian. Namely, 

Hi = - Y. O^- Y Oe (51) 

i,gy(£i) y(b) eg^(cd) 

where Ov,Oe are Pauli X operators with decorations of CZ operators on 2-links of v and 1-links of e respectively. 
The system has three types of Z 2 symmetry operators. Global 0-form symmetry operators are given by S'U) = 
n„gy(a) while a global 1-form operator is given by where is a 
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codimension-1 closed manifold intersecting with cd-edges. 

We study the boundary mode of the Hamiltonian. Consider a two-dimensional boundary which consists only of 


vertices of a,c, d (Fig. 15). The boundary can be seen as a three-colorable graph dK with color labels a,c,d where 


qubits live on cd-edges and a-vertices. The action of symmetry operators in the low-energy subspace can be written 
in terms of dressed boundary operators as follows: 


gia) 




Sib) 


eeM 


n 

(*d)6A 


(52) 


where is a closed loop on the boundary, Vi°'^ is the set of vertices of color a on the boundary and A is the set of 
all the 2-simplex on the boundary graph 9A. The following quantum critical Hamiltonian commutes with symmetry 
operators: 


Hboundary = - X ^ X, - ^ X, - 2^ A, [[ Z^. (53) 

y^V(a) e^£:(cd) y^yia) eGl-link(i;) iiGO-link(e) 

Note that this Hamiltonian is identical to a summation of trivial and non-trivial two-dimensional SPT Hamiltonians 
with 0-form and 1-form Z 2 ® 'Ll symmetry. This critical Hamiltonian can be transformed into two decoupled copies 
of critical Ising model by a duality transformation. 



FIG. 15: Boundary mode for a (0, 0, l)-form SPT phase. The left side of the boundary consists of vertices of color a, c, d while 
the right hand side consists of vertices of color a, b. A plaquette-like symmetry operator of control-Z operators is also shown. 


Next, let us consider a two-dimensional boundary which consists only of vertices of a, 6, c (Fig. 15). The boundary 
graph dK is a three-colorable graph with color labels a, 6, c where qubits live on a-vertices and 6-vertices while no 
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qubit lives on c-vertices. The actions of symmetry operators are 

•uGVt'*) t;GV<'’) 


n cz,,, 


(54) 


where is a plaquette-like product of CZ operators as shown in Fig. 15 The following quantum critical 

Hamiltonian, graphically shown, commutes with symmetry operators: 



With a bit of calculations, one can show that this critical Hamiltonian can be transformed into two decoupled copies 
of critical Ising model by a duality transformation. Thus, regardless of the choice of boundaries, one obtains two 
copies of critical Ising model or spontaneous breaking of Z 2 0 Z 2 symmetry. 

By gauging this three-dimensional SPT wavefunction in four dimensions, one can construct a non-trivial domain 
wall. The system consists of two copies of the (l,3)-toric code and one copy of the (2,2)-toric code with a domain 
wall. Let us denote anyonic excitations as follows: 
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^point ’ ^point ’ 
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'^loop 


(56) 


where "^^L6rar^e}: ^membrane}, {e/oop, Wioop} exhibit non-trivial braiding statistics respectively. The 

domain wall mixes the above excitations in an intriguing way. Upon crossing the domain wall, electric charges remain 
unchanged: 
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(57) 


while magnetic fluxes transform into composites of magnetic fluxes and superpositions of electric charges. Namely, 
one has 
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(58) 


/Q'\ ( 2 ^ 

where eioop) is a membrane-like object which consists of superpositions of and e/oop, and is characterized 

by a two-dimensional SPT wavefunction with 0-form and 1-form Z 2 0 X 2 symmetry while ^^point) ^ loop-like 

object which consists of superpositions of and i® characterized by a one-dimensional SPT wavefunction 

with 0-form Z 2 0 Z 2 symmetry. While the system consists of three copies of the toric code which do not belong to the 
same topological phase, a gapped domain wall which mixes three copies of the toric code can be constructed. This 
observation hints rather rich possibilities of gapped boundaries and domain walls in topological phases. 


E. Generic recipe 

Here we summarize the construction of q-form SPT Hamiltonians. Consider a d-dimensional simplicial lattice which 
is d-l-l-colorable with color labels oi,..., a^+i. Consider a partitioning of a positive integer d-l-1 with positive integers 
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Pf 


d+ 1 =Pl +P2 + • • ■ +Pm, l<Pl<P2< ■■■ <Pm (59) 

Let qj = pj — \ and Rj = pi + ... + pj-i. We will construct a model with (gi, < 72 , ■ • ■, 9 m)-form Z 2 ®"* symmetry. We 
place qubits on gr^-simplices of color kj := The trivial Hamiltonian is 

m 

Ho = -J 2 E (60) 

i=i 

where represents the set of all the gr^-simplexes of color kj. The non-trivial Hamiltonian Hi is constructed by 

applying the following finite-depth quantum circuit to the trivial Hamiltonian: 

C7(91.92.....9^) = Yl (61) 

where A represents the set of all the d-simplexes. The resulting Hamiltonian is 

m 

^1 = -E E ( 62 ) 

i=i vgA^'^P 


where 
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A Qj-ioTm symmetry operator is 


(63) 


v^‘‘P gM^'^P 


(64) 


where is a codimension-^j closed manifold intersecting with (jj-simplexes. 

When these models are constructed on a manifold with boundaries, non-trivial protected boundary modes can 
be supported. We do not have complete characterization of boundary modes. For the cases where m = 2, there 
are two possible choices of boundary vertices, and boundary modes are either the {qi,d— 1 — ( 72 )-toric code or the 
{q 2 ,d — 1 — ( 7 i)-toric code, which are equivalent to each other under Hadamard transformation. So, physics of the 
boundary mode does not depend on the choices of boundary vertices. For m > 2, boundaries can support protected 
gapless modes. Due to topological nature of the theory (ie diffeomorphism invariance), we expect that physics of the 
boundary mode does not depend crucially on choices of boundaries. However, we do not have independent argument 
for it. 

If one gauges the d-dimensional SPT wavefunctions in c? -|- 1 dimensions, one obtains a gapped domain wall in a 
d + 1-dimensional system which consists of the {pj,d + 1 — pj)-tonc code for j = 1,... ,to. The system possesses 
-dimensional electric charges, denoted by Cj, and (d — 1 — gj)-dimensional magnetic fluxes, denoted by rrij. Upon 
crossing the domain wall, electric charged remain unchanged: Cj —)■ Cj while magnetic fluxes get transformed as 
follows: 


171 j y nij s(ei,..., Cj—1, Cj-i-i,..., Cm) 


(65) 
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where s(ei,..., e^-i, Cj+i,..., Cm) is a (d — 1 — ( 7 j)-dimensional superposition of 
ei,..., Cj-i, Cj+i,..., Cm- Namely, s(ei,..., ej_i, e^+i,..., Cm) can be characterized by a 
(d — 1 — gj)-dimensional SPT phase with (gi,..., gj_i, gj+i,..., gm)-form symmetry, 

magnetic fluxes and fluctuating charges will exhibit non-trivial multi-brane braiding statistics. 


electric charges 
wavefunction of 
We expect that 


IV. FAULT-TOLERANT LOGICAL GATE 

In this section, we comment on applications of our construction to the problem of classifying fault-tolerantly 
implementable logical gates in topological quantum codes. By topological quantum codes, we mean quantum error- 
correcting codes, supported on lattices, that can be characterized by geometrically local generators. Namely, we shall 
argue that all the examples of SPT phases with generalized global symmetry proposed in this paper have corresponding 
fault-tolerantly implementable logical gates in topological quantum codes living in one more dimensions. 


A. Logical gate and domain wall 

We begin by recalling the connection between classifications of fault-tolerant logical gates and gapped domain 
walls |28l I29j . The underlying difficulty in quantum information science is the fact that quantum entanglement decays 
easily and qubits need to be protected from noise and errors. In theory, this challenge can be resolved by using 
quantum error-correcting codes where single qubit information is encoded in many-body entangled states such that 
local errors do not destroy the original information. Then, one can perform quantum computation fault-tolerantly 
inside a protected subspace (codeword space) of a quantum error-correcting code by performing error-correction 
frequently. 

A naturally arising question concerns how to implement logical gate operations inside the codeword space. Ideally, 
one hopes to perform logical gates by transversal implementations of unitary operators which have tensor product 
form, acting on each qubit individually. For such a transversal logical gate, local errors do not propagate to other 
qubits, and thus its implementation is fault-tolerant. One may also fault-tolerantly implement logical gates which can 
be expressed as finite-depth local quantum circuits. However, if a logical gate implementation requires a highly non¬ 
local and complicated quantum circuit, local errors may propagate to the entire system in a uncontrolled manner. As 
such, it is important to find/classify fault-tolerantly implementable logical gates in quantum error-correcting codes 
To gain some intuition on the restriction on fault-tolerant logical gates, consider the two-dimensional toric code. 
The system has string-like Pauli X and Pauli Z logical operators which act non-trivially inside the ground state space 
(codeword space) and have transversal form. However, the toric code does not admit any other transversal logical 
gates except for products of string-like Pauli operators. In this sense, the toric code has a rather restricted set of 
transversal logical gates. Thus, we hope to find quantum codes with a larger set of fault-tolerant logical gates 

Let us look at an example where one can implement non-Pauli logical operators transversally. Consider a two- 
dimensional system which consists of two decoupled copies of the toric code. Namely, we assume that the first copy 
lives on a square lattice and the second copy lives on a dual lattice as shown in Fig. [T^a). The Hamiltonian is given 


^ There are logical gates which do not admit finite-depth circuit implementation, but can be implemented in a rather simple manner. For 
instance, a Hadamard-like logical gate can be implemented in the two-dimensional toric code by shifting the lattice sites in a diagonal 
direction, followed by transversal application of Hadamard operators. 

® However, a larger set of fault-tolerant logical gates often implies weaker error tolerance. See m 
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by 


H = - Y. 41 )- ^ 42 ) _ ^ 5 ( 1 )_ ^ ^( 2 ) (gg) 

„gv(i) jjGV'f^) pGP(i) pGP(2) 

where two qubits are placed at each site. Here represent the sets of vertices on lattice 1 and 2 respectively 

while represent the sets of plaquettes. In this system, one can apply a logical control-Z gate among two copies 

of the toric code transversally. Specifically, we apply control-.Z gates to pairs of qubits in the first and second copy at 
the same sites. This transversal operation preserves the ground state space, but has non-trivial action on it. Recalling 
the transformation rules under conjugations by control-.Z gates, CZ{Xi)CZ = (Xi)Z 2 , CZ(X 2 )CZ — Zi(X 2 ), one 
notices that transversal control-Z gates transform anyonic excitations as follows: 


Cl —>■ Cl mi —>■ e2mi 62 —>■ 62 m2 —>■ eim2. 


(67) 


In general, membrane-like fault-tolerant logical operators transpose labels of anyonic excitations in two-dimensional 
topological quantum codes [ 55 ] . 


(a) lattice 1 lattice 2 
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FIG. 16: (a) Two copies of the toric code. The first copy lives on a square lattice with solid lines while the second copy lives on 
a dual lattice with dotted lines. Grey dots represent two qubits from each copy of the toric code. Control-gates are applied 
to pairs of qubits at the same sites, (b) A gapped domain wall constructed from a control-Z logical gate. We apply the logical 
gate only on the right hand side of the lattice. 

In order to construct a gapped domain wall, let us split the entire system into the left and right parts and apply 
the transversal control-Z gate only on the right hand side of the lattice. This transforms the Hamiltonian into the 
following form: 


P — Hleft P Pwall P Pright 


( 68 ) 


where Pieft and Pright remain unchanged while Pwaii can be viewed as a gapped domain wall which connects 


Pleft and Pright- Upon crossing the domain wall, anyonic excitations are transposed according to Eq. ( [67| ). As 
this observation implies, given a d-dimensional non-trivial fault-tolerant logical gates in a d-dimensional topological 
quantum code, one can construct a corresponding domain wall since non-trivial logical gates would transform types 
of excitations. 
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B. Logical gates and generalized global symmetry 

Next, let us establish the connection between SPT phases and gapped domain walls in the context of fault-tolerant 
logical gates. For this purpose, we revisit the gapped domain wall in the two-dimensional toric code. Consider 

a two-dimensional trivial symmetric system where a one-dimensional SPT wavefunction with Z 2 ® Z 2 symmetry is 
inserted as in Fig. ©a). By gauging the entire system with respect to symmetry, one obtains two copies of the 

toric code with a gapped domain wall as shown in Fig. |17[ b). Let 7 be a one-dimensional line where a one-dimensional 
SPT wavefunction is placed. Observe that one can move one-dimensional SPT wavefunction to a different location 7 ' 
by applying a symmetric local unitary transformation U which acts only on qubits enclosed by 7 and 7 '. By gauging 
the system, this process of moving one-dimensional SPT wavefunction is equivalent to moving a domain wall from 7 
to 7 ' by applying a local unitary transformation U on qubits enclosed by 7 and 7 '. By sweeping the domain wall over 
the entire system, one can implement a non-trivial logical gate. Namely, this implements the control-^ gate among 
two copies of the toric code. In this sense, a one-dimensional SPT wavefunction with Z 2 ® Z 2 symmetry corresponds 
to the control-Z logical gate in two copies of the two-dimensional toric code. Here we would like to emphasize that 
the logical gate acts on a system with intrinsic topological order (the Z 2 ® Z 2 toric code) while the logical gate was 
constructed by gauging the one-dimensional Z 2 ® Z 2 SPT wavefunction. 

In |28] , we employed this idea to construct fault-tolerantly implementable logical gates in the d-dimensional quantum 
double model. Namely, for the d-dimensional quantum double model with hnite group G, the domain wall can be 
constructed by gauging a {d — l)-dimensional SPT wavefunction in d dimensions. In other words, one is able to 
construct logical gates, gapped boundaries and domain walls by using d-cocyle functions in the d-dimensional quantum 
double model. We then demonstrated that, for non-trivial domain walls, there exist corresponding non-trivial logical 
gates which can be implemented by hnite-depth local quantum circuits. It was also found that the non-triviality of 
domain walls can be verihed by computing slant products sequentially. 

(a) (b) 




FIG. 17: A fault-tolerant logical gate constructed from SPT wavefunctions and gapped domain walls, (a) Moving a one¬ 
dimensional SPT wavefunction by finite-depth symmetric quantum circuits, (b) Moving a gapped domain wall by finite-depth 
quantum circuits. 

A similar argument applies to SPT phases with generalized global symmetry. Consider a d-dimensional SPT phase 
with (gi,..., (Z 2 )®’" symmetry. By gauging this SPT phase in d-I- I dimensions, one obtains a gapped 

domain wall in a d -I- 1-dimensional topological phase. This d -I- 1-dimensional system consists of {pj, d + 1 — pj )-toric 
code for j = 1,..., m. One can move (d — l)-dimensional SPT wavefunctions by applying symmetric finite-depth 
quantum circuits. This implies that one can sweep the domain wall over the entire system by applying a finite- 
depth local unitary circuit. With some speculation, one notices that this implements a logical gate among 
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{pj,d+ 1 — pj)-toric code. For instance, by using a three-dimensional SPT phase with (0, 0, l)-form symmetry, one 
can construct a G®'^Z logical gate acting among two copies of the (l,3)-toric code and one copy of the (2,2)-toric 
code in four dimensions. It is interesting to observe that logical gates can be implemented among several copies of 
the toric code which belong to different quantum phases. Our findings hint rich possibilities of fault-tolerant logical 
gates, as well as gapped domain walls, in topological phases of matter. 


V. DISCUSSIONS 

We will conclude the paper with discussions on generalizations of the model and implications of our results. 


A. Generalizations 

While our treatment has been limited to systems with Z 2 symmetries, the construction can be generalized to systems 
with arbitrary abelian symmetries. Here we illustrate the idea for Zjv symmetry. We define the following generalized 
TO-qudit control-Z gate 

{/(m) = ^ exp j | 5 i,...,g„)( 5 i,...,g™| (69) 

where gj = 0,..., N — 1. For N = 2, this reduces to the multi-qubit control-Z gate. To construct SPT phases with 
0-form symmetries, consider a (d-l- 1) colorable graph A in d dimensions and assign qudits (Wstate spins) to vertices. 
On a colorable graph, one is able to assign parity F’(A) = ±1 to each d-simplex such that neighboring simplexes have 
opposite parity signs |33] . We shall apply the following local unitary to the trivial symmetric Hamiltonian 

U (70) 

A 

where acts on d-I-1 qudits on the d-simplex A. The system possesses d-|-1 copies of Z^r symmetries, associated 

with each different color label ai,... , 0 ^+ 1 . To construct SPT phases with higher-form Zjv symmetries, one places 
qudits according to the partition of color labels. 

Our construction of SPT phases does not exhaust all the possible bosonic SPT phases. Yet, by changing the choices 
of symmetry operators, one is able to construct some other SPT phases. Let us illustrate the idea by considering two- 
dimensional SPT phases with Z 2 symmetry. The proposed model possesses three copies of Z 2 symmetries, captured 
by three symmetry operators S'a 0 S'b ® Sc, associated with three color labels A, B, C. It is possible to view the model 
as an SPT Hamiltonian with one copy of Z 2 symmetry by imposing SaSbSc as the single Z 2 symmetry operator. 
This reduces the model to the one proposed by Levin and Gu m- One may also choose to impose two copies of Z 2 
symmetries by using Sa O SbSc- Thus, there are three possible SPT phases associated with symmetries: 

SaSbSc Sa ® SbSc Sa'^ Sb Sc- (71) 

Recall that, in two spatial dimensions, SPT phases protected by abelian symmetries can be classified into three classes, 
call type-I, type-H and type-HI [2]. Levin and Gu showed that the SaSbSc model corresponds to the type-I model. 
In [ 21 ], we showed that the Sa Z> Sb G> Sc model corresponds to the type-HI model by the domain wall argument. It 
is an interesting question to determine the type of the Sa 0 SbSc model. 
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B. Previously known models 

We comment on the relations between our models and previously known models in three dimensions. Our model 
can be characterized by partitions of the integer 4 (which is the space-time dimension). Namely, possible models 
are characterized by (1,1,1,1), (2,1,1), (2,2) and (3,1). The (1,1,1,1) model possesses 0-form SPT order with 
S'a ® S'b 0 S'c' O S'i). One can pick different choices of symmetries, such as SaSb O S'c ® 5'd, to construct different 
classes of 0-form SPT phases. The (2,2) model possesses 1-form symmetries only, and seems to be described by 
a theory containing BAB term after gauging. Such TQFTs typically belong to a (rather simple) subclass of the 
Walker-Wang model m- The (3,1) model involves 2-form symmetries which have not been discussed much in the 
literature. However, we think that the model can be reduced to a known one by a duality transformation which 
exchanges charges and fluxes (in other words, Pauli X and Z operators). The (2,1,1) model has one copy of 1-form 
symmetry and two copies of 0-form symmetries, denoted by Sab Sc 'Z)Sd- Consider a realization with Sab SqSd ■ 
After gauging, the model seems to be identical to the Birmingham-Rakowski model with Z 2 symmetry introduced 
in [46] . The Birmingham-Rakowski model can be discussed in a more generic framework, called the Mackaay TQFT 
where degrees of freedom are placed both on edges and faces as in our construction [47]. As for the model with 
Sab Sc Sd symmetry, we were not able to find similar constructions in the literature. Whether this construction 
is truly new or not waits further verifications. 

Kapustin and Thorngren utilized the notion of 2-group to construct SPT phases protected by 0-form and 1-form 
symmetries on lattices. An important difference between our approach and the 2-group construction is that the 
2-group construction possesses non-fiat connections while our models consist only of fiat connections. In simpler 
words, our models modify X-type vertex terms only while the 2-group construction modifies Z-type plaquette terms 
too. However, duality transformations often allow one to construct models with non-fiat connections from the ones 
with flat connections. Indeed, some of the 2-group constructions are equivalent to models only with flat connections. 
Similarly, some of our models are equivalent to models with non-flat connections via duality transformations. Such 
complementary viewpoints may allow one to construct further examples of interesting TQFT models. At this moment, 
we were not able to fit our (2,1,1) model into the 2-group construction despite the fact that the model has 0-form and 
1-form symmetries only. In |21| , Kapustin and Thorngren briefly mention possible generalizations using the notion of 
q-group (g > 2). Presumably, such constructions will involve several symmetry operators of different dimensionality, 
and our construction may give concrete examples of such generalizations. 


C. Future problems 

Other questions and future problems are listed below. We did not discuss physical properties of gauged models in 
depth. Upon gauging SPT Hamiltonians, one typically obtains topologically ordered Hamiltonians whose braiding 
statistics are twisted due to the decorations added to matter fields. We expect that gauged versions of our models 
will exhibit rather exotic topological order which may be beyond known theoretical frameworks. For one thing, 
in two dimensions, the 0-form Z 2 0 Z 2 (8) Z 2 SPT phase considered in this paper has been shown to be dual to a 
non-abelian topological phase upon gauging HIMIIISI- This hints a possibility of interesting non-abelian statistics 
which involve both particles and loops by gauging higher-form SPT phases. SPT phases with g-form global symmetry 
provide a number of interesting quantum critical Hamiltonians as boundary modes. Analytical and numerical studies 
of such boundary modes may provide further insights into problems of quantum criticality in higher dimensions. 
Spatial dimension of symmetry operators can be non-integer [451 149] . Namely, one can construct an SPT Hamiltonian 
protected by fractal-like symmetry operators. Studies of such fractal SPT phases and their gauged models may be an 
interesting future problem with applications to efficient magic state distillations. One drawback of our approach is that 
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the proposed models do not have full diffeomorphism invariance due to the use of colorable graphs. While it is possible 
to coarse-grain or fine-grain graphs by retaining colorability |34j . the full verification of fixed-point properties is an 
important future problem. Finding field theoretical descriptions of proposed model is also an interesting project Ell- 
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